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We provide a simple scheme for predicting the electronic structure of van-der-Waals bound sys- 
tems, based on the mere knowledge of the electronic structure of the subunits. We demonstrate this 
with the example of nano-peapods, consisting of polythiophene encapsulated in single- wall carbon 
nanotubes. Using density functional theory we disentangle the contributions to the level alignment. 
The main contribution is shown to be given by the ionization potential of the polymer inside the 
host, which, in turn is determined by the curvature of the tube. Only a small correction arises 
from charge redistributions within the domains of the constituents. Polarization effects turn out to 
be minor due to the cylindrical geometry of the peapods and their dielectric characteristics. Our 
findings open a perspective towards designing opto-electronic properties of such complex materials. 
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The creation of nano-structures with tailored proper- 
ties is an emerging field, opening exciting perspectives in 
all kinds of research areas and technologies. Exploiting 
the potential of nano-technology requires a thorough un- 
derstanding of the respective building blocks as well as 
their interplay. Only that way, specific novel features 
can be designed. Theory thereby plays an important 
role in getting deeper insight. For a quantitative anal- 
ysis one needs to make use of the most sophisticated ab- 
initio methods, which are in most cases too involved to 
be applied to real systems. Hence, simple schemes to 
describe specific properties of a nano-structure from the 
knowledge of its constituents, are most valuable tools to- 
wards materials design. In this work, we propose such a 
scheme to predict and hence tailor the electronic struc- 
ture of weakly bound nano-hybrids. We demonstrate our 
approach with peapods - carbon nanotubes {pods) accom- 
modating molecules or atoms (peas) in their cavities. 

Nano-peapods have attracted scientific interest since 
their discovery.!^ They combine the unique mechanical 
properties of the nanotube with the opto-electronic prop- 
erties of the guest. This way, very stable hybrid materi- 
als with tunable characteristics can be formed. In most 
applications, the role of the guest is to alter the elec- 
tronic structure of its host in order to achieve the de- 
sired features. For example, they have been proposed for 
opto-electronics,!^ electron-transport devices,^ and even 
for more fundamental purposes, like one-dimensional spin 
arrays with an outlook towards quantum computing.^ 

While nanotubes with fullerenes inside have been stud- 
ied extensively over many years, organic molecules in- 
side the cage represent a more recent topic. After the 
encapsulation of Zn-diphenylporphyrin,l^ tetracyano-p- 
quinodimethane,^ perylene derivates,'^ and /3-carotene,^ 
the successful synthesis of nano-peapods with thiophene 
oligomeres was reported.El The particular feature of the 
latter is their ability to emit light in the visible range of 



the spectrum. This provides new opportunities for opto- 
electronic devices as pure carbon nanotubes are infra-red 
emitters, while being passive in the visible range. The 
question immediately arising is how the electronic and 
hence also the optical properties of such systems can be 
designed and predicted prior to synthesis. 

To address this point we exemplarily investigate a se- 
ries of peapods, consisting of (n,0) single-wall carbon 
nanotubes (SWNTs) and polythiophene representing the 
guest. In the following, we use the short-hand notation 
pT@(n,0), and (oo,0) refers to a single graphene sheet 
as the limiting case of a tube with infinite diameter. We 
point out, that these systems are purely van der Waals 
(vdW) bound^ ^ as will be seen in more detail below. 

All calculations presented here are based on den- 
sity functional theory (DFT) using ultra-soft pseudo- 
potentials as implemented in the program package 
QUANTUM-ESPRESSO.io A plane-wave cutoff of 30 
Ry is chosen, and a total of 4 k points is taken into 
account for the one-dimensional Brillouin zone integra- 
tions. We adopt a super-cell approach with a vacuum 
size of 8 A to avoid interactions between adjacent nan- 
otubes. Exchange-correlation effects are treated within 
the generalized gradient approximation (GGA) in the 
PBE parameterization.^^ For geometry relaxations, the 
vdW density functional (vdW-DF)^^ is employed in a 
post-scf manner. A self-consistent vdW-DF procedure- 
for selected systems proved the charge density not to 
be affected, as already found by Thonhauser for other 
systems. To ensure commensurate lattice parameters 
along the tube axis, we double the repeat unit of the zig- 
zag SWNT and stretch the polymer by about 8%. We 
consider two different geometric arrangements, one with 
the polymer centered inside the tube and another one 
where it is sitting at its optimal distance from the tube 
wall. 

In Fig. [l] the band structures obtained by DFT cal- 
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FIG. 1. Band structures of the peapods pT@(n,0) with n 
= (12, 13, 14, 22) shown in the vicinity of the Fermi energy 
(black horizontal lines). Full turquoise lines, thick gray lines, 
and green dots represent the electronic bands of the isolated 
polythiophene, isolated nanotube, and the peapod, respec- 
tively. All band energies are aligned with respect to the vac- 
uum level. The peapod bands originating from the polymer 
appear to be rigidly shifted by an amount AIP compared to 
the isolated polymer, as indicated by the arrows. 



culations of the isolated polythiophene (turquoise lines), 
the isolated nanotubes (bold gray lines), and the corre- 
sponding peapods (green dots) are plotted for four zig- 
zag SWNTs. For any system consisting of fully non- 
interacting subunits, the total electronic structure is a 
superposition of those of its constituents, aligned with 
respect to the vacuum level, i.e., the electrostatic po- 
tential at infinite distance. This so-called Schottky-Mott 
limit^^ should generally apply to vdW bound systems, 
where no charge transfer between the constituents oc- 
curs. Hence, the bands in Fig. [l] are aligned in this man- 
ner. In fact, the electronic states of both subunits re- 
tain all their features when forming the band structure 
of the peapod. We observe, however, that those pea- 
pod bands which are originating from polymer states ap- 
pear rigidly shifted downward (by an energy labeled AIP) 
compared to the band structure of the isolated polythio- 
phene. These shifts are rather small, ie., 0.12, 0.14, 
0.15, and 0.08 eV for pT@(12,0), pT@(13,0), pT@(14,0), 
and pT@(22,0), respectively. In the following, we aim 
at explaining their origin, thereby enabling a prediction 
of the peapod's band structure from the knowledge of 
its constituents' bands alone. For this, we have to con- 
sider three different factors. The first one is specific to 
carbon nanotubes, the second one emerges from charge 
rearrangements between weakly interacting entities, and 
the third factor is connected with mutual polarization of 
the subunits. These will be discussed in the following. 
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FIG. 2. (a) Charge density of a (16,0) SWNT as a function of 
the radial coordinate r, averaged along the tube axis and the 
azimuthal angle. The blue and red areas below the curve rep- 
resent the charges inside and outside the tube perimeter, Nin 
and Nout, respectively, (b) The corresponding electrostatic 
potential averaged along the tube axis, (c) Difference in the 
electrostatic potential inside and outside the tube perimeter 
as a function of the normalized charge asymmetry. 



The first factor stems from the curvature of the car- 
bon nanotubes. It gives rise to a charge asymmetry be- 
tween the inside and outside area with respect to the 
tube wall. This is visualized in the top panel of Fig. [2^. 
Here, the charge density, averaged along the tube axis 
and the azimuthal angle, is depicted for a (16,0) SWNT 
as a function of the distance from the tube center, r. The 
areas below this curve left and right of the tube perime- 
ter (gray dashed line), marked in blue and red, respec- 
tively, represent the corresponding number of electrons. 
The observed charge asymmetry results in an offset in 
the electrostatic potential inside the host compared to 
the vacuum level. This is illustrated in Fig. ^ where 
we plot the potential difference = ^out — ^in foi" a 
series of (n,0) tubes versus {Nout — Nin)/Ntot' We find 
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to depend linearly on this normalized charge asym- 
metry and to be proportional to the tube curvature Xjd. 
As a consequence of this electrostatic effect, molecules 
positioned inside or outside the tube experience a differ- 
ent electrostatic potential, and thus the molecular levels 
move downwards as compared to the isolated molecule. 

The second factor originates from a small charge re- 
distribution when the guest is getting in contact with the 
host. To quantify its impact, we evaluate the difference 
density Ap = ppeapod - Ptube - Ppoiymer, whcrc the densi- 
ties of the polymer, ppoiymen and that of the tube, ptuhei 
have been computed in the same geometry as that of 
the peapod, ppeapod- This charge density difference is de- 
picted in the upper part of Fig.jsjfor pT@(18,0) (left) and 
pT@(14,0) (right) for the polymer in the wall and center 
position, respectively. We analyze the charge transfer be- 
tween the two subunits by integrating Ap separately over 
the regions of the polymer and nanotube, respectively, 
which we define through the condition V ppeapod = as 
indicated by the green ovals. As expected, the so ob- 
tained charge, AQ, exhibits a systematic decrease as the 
tube size increases, reaching zero for the centrally situ- 
ated polymer, while it approaches the tiny value of 0.004 
transferred electrons in case of the optimally spaced sub- 
systems. Hence, we can conclude that there is negligible 
charge passing from one subunit to the other as expected 
for vdW bound systems. The main charge rearrangement 
takes place within the individual domains leading to for- 
mation of electric dipoles as seen in Fig. [3| To explore 
the impact of such charge dipoles quantitatively, we solve 
Poisson's equation V'^^corr = Ap for the differential den- 
sity Ap. This yields another contribution, $corr, to the 
electrostatic potential which is analogous to the bond 
dipole effect discussed in the context of molecule/metal 
junctions. 

The results obtained so far are summarized in Fig. |4] 
Here, the ionization potential (IP) of the polymer inside 
the SWNT as a function of the tube diameter is displayed 
for the two different positions of the pea inside the pod. 
The turquoise lines correspond to the IP of the isolated 
polymer, and the green filled circles represent the DFT 
results for the peapod. Their difference is identical to the 
quantity AIP indicated by the arrows in Fig. [l] It pro- 
vides an estimate for the accuracy we can expect when 
considering the Schottky-Mott limit as a zeroth approxi- 
mation to the level alignment. The two factors discussed 
above represent corrections to this picture. First, the 
fact that the polymer experiences a different IP inside 
the pod compared to outside, gives rise to the correc- 
tion A<l>. The resulting data are marked by gray crosses. 
Including this correction, determined by the tube curva- 
ture, leads to very good agreement with the results for 
the combined system for realistic tube diameters as ac- 
complished experimentally.^ The deviation is within the 
range of O.leV. This is an exciting result as only the 
knowledge of the pristine systems is required to obtain 
this approximate IP. 

Evaluating these approximate IPs for the central and 
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FIG. 3. Top: Charge density difference integrated along the 
tube axis with isovalues of +0.0001 (red) and -0.0001 (blue) 
electrons per a.u.^. Results are displayed for pT@(18,0) (two 
pictures at the left) and pT(14,0) (two pictures at the right) 
for central and optimal position of the polymer, respectively. 
The green ovals mark Vp = 0, separating the domains of 
the two subsystems. Bottom: Charge exchanged between the 
nanotube and the polymer in dependence of the tube diam- 
eter. The filled spheres and the triangles indicate the cases 
of polythiophene in the center of the tube and its optimal 
distance from the tube wall, respectively. 



the off-centered position further, the deviations from the 
peapod data are somewhat larger for the latter. This can 
be understood in terms of guest-host interaction, which is 
also responsible for the pronounced discrepancies found 
for very small tubes. Here, the corrections arising from 
the charge dipoles, ^com come into play. Taking them 
into account (orange triangles), leads to excellent agree- 
ment with the DFT results for both types of positions 
and the entire curvature range. Obviously, for very nar- 
row tubes, the charge rearrangement upon encapsulation 
is dominant over the curvature effect. One should point 
out, however, that the formation of such peapods is en- 
ergetically unfavorable.'^ 

At this point, we should assess the reliability of DFT 
in describing such scenarios as discussed here and ex- 
ploiting its predictive power. To start with, we point 
out again that we are dealing with vdW-bound systems. 
It is clear that conventional (semi) local functionals can 
fail badly in this context. However, vdW- based 
on the PBE density has proven successful for the sys- 
tems under investigation^ without leading to spurious 
charge transfer. Concerning the DFT band-gap problem, 
we point out that for our analysis we do not make use 
of band gaps. We only evaluate the valence band struc- 
tures to demonstrate the principles behind the alignment 
of the respective HOMO levels. The suggested procedure 
can equally be based on calculations performed with hy- 
brid functionals or the Gl^ approach, or even solely on 
experimental data. 
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FIG. 4. Ionization potential as a function of the inverse 
tube diameter for the polymer in the center of the tube (top) 
and its optimized position (bottom) . The green filled symbols 
indicate the IP obtained by DFT calculations for the peapod. 
The turquoise horizontal lines depict the IP of the isolated 
polymer. Gray crosses show the IP of the polymer inside the 
cavity at its respective position. The orange open symbols 
mark the corresponding values when accounting for the charge 
dipoles arising between the polymer and the tube. 



Evaluating the third factor, however, which is con- 
nected with the polarization response of the SWNT 
upon adding a hole (electron) in the polymer's HOMO 
(LUMO), (semi)local DFT would certainly fail. We em- 
phasize, however, that this effect is small and, can equally 
be understood in terms of electrostatic considerations. 
Such scenario has been discussed for the example of ben- 
zene physisorbed on graphite, showing that the GW 
band gap of benzene is substantially reduced upon ad- 
sorption, which could be well interpreted by a simple 
model based on an electrostatic image potential. To 
quantitatively estimate these effects for our system, we 
adopt two different electrostatic models representing two 
limiting cases. 

The first model^^ considers a point charge at the equi- 
librium polymer-SWNT distance of around 3.8 A in front 
of a dielectric slab of thickness t. The results are shown in 
Fig. [5^, where the electron energy E of the point charge 
with the induced image charge is plotted as a function of 
t for various dielectric constants s of the slab material. 
Approximating the thickness of graphene by its electron 
density's full-width-at-half-maximum value, represented 
as vertical dashed line, we observe a strong dependence of 
the image-potential effect on the dielectric properties of 
the slab. A perfect metal {e = oo) would result in level 
shifts of almost 1 eV - independent of the slab thick- 
ness. This corresponds well to GW results for the fron- 
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FIG. 5. (a) Image charge model of a dielectric slab in the 
presence of an elementary charge: Electronic energy shift as a 
function of slab thickness t for various dielectric constants s. 
The vertical dashed line at tgr indicates the effective thickness 
of a graphene sheet, (b) Image charge model of an infinite 
metallic cylinder hosting an elementary charge placed at a 
distance of 3.8 A from the cylinder wall: Electronic energy 
shift as a function of the inverse cylinder radius 1/R. In both 
graphs the same color code for the dielectric constants is used. 



tier molecular orbitals of benzene adsorbed on graphite.!^ 
(Note that in these calculations, the adsorption distance 
was much smaller, 3.25 A.) However, the SWNTs stud- 
ied in this work are less polarizable (e.g. 5 for (14,0) 
tube), and, therefore, polarization shifts of about 0.1-0.3 
eV only are expected (see Fig. [5|i). 

The above estimates are based on a planar dielectric 
slab. In order to take the curvature of the SWNT into ac- 
count, we employ a second model where a point charge is 
located inside an infinitely long metallic cylinder, and 
evaluate the energy £^ as a function of inverse cylinder 
radius. As follows from Fig.[5]3, metallic cylinders of radii 
corresponding to (22,0), (16,0), and (10,0) SWNTs would 
give rise to polarization shifts of -0.4 eV, -0.2 eV, and 
eV, respectively (gray line). Note that for (10,0) the 
charge resides in the center of the tube where no polariza- 
tion shift is to be expected to due symmetry reasons. As, 
however, the limit of metallic cylinders present an upper 
bound to the polarization shift, a considerable reduction 
is to be expected, by considering the finite dielectric con- 
stant of the SWNT. Unfortunately, no analytical solution 
is available for a point charge off the axis inside cylindri- 
cal dielectric of finite thickness. In order to estimate the 
reduction to the finite e starting from the exact results 
for a perfect metallic cylinder, we therefore proceed as 
follows. We adopt the polarization shifts from the pre- 
vious model (point charge in front of a planar dielectric 
slab) and include these values also in Fig. [sJd for zero 
curvature {1/R = 0), i.e., the planar case. We further 
assume that the curvature dependence for finite dielectric 
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constant exhibits the same behavior as that for the per- 
fect metal. With this extrapolation, we can now estimate 
the image charge effect for arbitrary nanotube diameter 
and given dielectric constant as presented by the dashed 
curves in Fig. [sJd. This reasoning allows us to estimate 
the overall polarization effects to be in the order of 0.1 
eV or less for SWNT diameters and values relevant for 
our examination here. 

To summarize, we have discussed the level alignment 
of weakly bound nano-peapods. We find that due to 
the vdW bonding, there is negligible charge transfer be- 
tween the subunits, preserving the nature of the single 
bands within the electronic structure of the combined 
system. The alignment of the HOMO levels is governed 
by the ionization potentials of the separated entities with 
a small correction arising from charge rearrangements in- 
duced by their proximity. It is important to consider the 
IP of the pea inside the tube, accounting for the dif- 
ferent electrostatic environments inside and outside its 
perimeter. These IPs are uniquely determined by the 



tube curvature. We provide them for a series of SWNTs 
in this work and encourage their measurement, e.g.^ by 
core level or photoemission spectroscopy. Polarization 
could be an important factor und needs, in general, to 
be taken into account. However, due to the cylindrical 
symmetry of nano-peapods and their dielectric proper- 
ties, the corresponding electronic energy shifts are con- 
siderably reduced. 

Our findings allow for predicting the band structure 
of nano-hybrids from the mere knowledge of its con- 
stituents. From our examples, which are representative 
for molecular systems used in opto-electronics, we esti- 
mate the error bar to be of the order of 0.1 eV. Such pre- 
dictions could be extremely useful for designing electronic 
properties of nano-structures towards opto-electronic de- 
vices or nano-transistors where tuning the electronic lev- 
els is crucially important. 

Work carried out within the NanoSci-ERA project 
NaPhoD (Nano-hybrids for Photonic Devices), supported 
by the Austrian Science Fund, projects 1107 and S9714. 
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